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ABSTRACT 

Forthcoming datasets from the Planck experiment and others are in a position to probe the CMB non- 
Gaussianity with higher accuracy than has yet been possible and potentially open a new window into the 
physics of the very early universe. However, a signal need not necessarily be inflationary in origin, and 
possible contaminants should be examined in detail. One such is provided by early universe magnetic 
fields, which can be produced by a variety of models including during an inflationary phase, at phase 
transitions, or seeded by cosmic defects. Should such fields have been extent in the early universe they 
provide a natural source of CMB non-Gaussianity. 

Knowledge of the CMB angular bispectrum requires the complete Fourier-space (or "intrinsic") bis- 
pectrum. In this paper I consider in detail the intrinsic bispectra of an early-universe magnetic field for a 
range of power-law magnetic spectra. 



1. Introduction 



Observations of the cosmic microwave background (CMB) strongly suggest that the dominant source of pertur- 
bations in the early universe was adiabatic and nearly Gaussian in nature ( Komatsu et al.|[20T0 l, consistent with an 
inflationary scenario. Inflationary scenarios also induce a small amount of non-Gaussianity, typically characterised by 
a parameter /nl ( [Liguori et al.|2010] for example) which from the WMAP seven-year results is currently constrained to 



be - 10 < /nl < 74 at 95% confidence ( Komatsu et al. 2010 1. It is expected that these bounds will tighten significantly 
with the forthcoming Planck data ("Liguori et al. 2010 1. A detection of a significant /nl could provide vital information 
about the early universe and high-energy physics, as it could rule out a variety of inflationary models including the 
simplest single-field models. 

However, the inflationary model does not exclude the possibility that non-linear sources might also play a role 
in sourcing perturbations. Frequently-studied examples of such sources include cosmic defects ( Turok et al.|1998 for 
example) and cosmological magnetic fields ( Giovannini|2004 Tsagas|2005 for example). 



Magnetic fields are observed on many scales in the cosmos, including on cluster scales with a coherence length on 
the order of megaparsecs and field strengths on the order of between nano-Gauss and micro-Gauss (I Kronberg|1994 



Grasso & Rubinstein||200l[ Giovannini||2004[ Xu et al.||2006[ Kahniashvili et al.]|2010| . Fields on larger scales are 
notoriously difficult to detect but there are some suggestions that such fields might exist with field strengths up to 
the order of micro-Gauss. Recent observations of radiation from blazars implies a lower limit in the extra-galactic 
medium of Z? > (9(10"'^)-(9(10"'^), permeating extracluster voids (???). The presence of fields on such fields implies 
either a primordial origin or an efficient transfer of fields from within galaxies deep into the intergalactic medium. The 
precise origin of fields on such large scales remains debated but magnetogenesis scenarios exist which were viable 
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in the extremely early universe. Such scenarios can be produced directly during inflation ( Turner & Widrow 1988J 
[Bassett et ai:][250T] IProkopec & Puchwein"2004'; "Giovannini '2008 ; Campanefli et al.' 2008'; Bamba et al.1|2008i ) or 
at a phase transition ( ,Baym et al. 1996^ Martin & Davis 1995; Hindmarsh & Everett^ 1998; Boyanovsky et al.|2003 



Kahniashvili et al.|2009b i. Fields can also be produced by the production of non-linear vorticity from linear density 



perturbations( |Gopal & Seflii|2005t 'MataiTes eet al.|2005[ [Takahashi etal.|2005[ [Ichiki et al.|2"006l [Siegel & Fryl^OOe 



Kobayashi et al.|2007[ Maeda et al.||2009j ), but the impact these have on the CMB is complicated by their evolving, 



non-trivial nature. 



The impact primordial magnetic fields have on the CMB and its anisotropics have been well-studied, with (Bar- 
|row et al.|1997l|Subramaman & Barrow] 1998bl|Dun-er et al.|1998l|Koh & Lee|2000l [Kahniashvih et al.|2000l [Mack| 
let al.|2002trClarkson et al.|20031[Lewi7 2004; K ahniashvih"&R atra 2007, Kahniashvil T et al.|2009al |Yamazak retaL] 
[2008','Fineni et al.'2008'' Paoletti et al.|2009nBonvin & Caprini|2010,|Giovannini,2010(|^mazaki et al.|2010nPa'oietti 



[& Finelli 2010; Kahniashvih et al.|2010 1 being some instructive examples. While older literature tended to assume a 
homogeneous background component with an inhomogeneous perturbation, more recent work has typically focused 
on tangled configurations without a background component and I assume this throughout. These studies fairly consis- 
tently suggest that the field is constrained to be of at most nano-Gauss in magnitude. The spectral index is restricted 
to be approximately scale-invariant ( Yamazaki et al.|[2010[ Paoletti & Fiiiem||2010| l, with limits growing extremely 
tight for a primordial magnetic field with index far from scale-invariance ( Caprini & Durrer]|2001[ l. A large-scale 
homogeneous field also introduces characteristic correlations between multipole moments with Al e {-2, 0, 2) and 
Am e {0, +1, ±2) which vanish in the standard scenario ( [Kahniashvili et al.,2008j . 

However, the magnetic 2-point signal is overwhelmed on large-scales by the standard perturbations, with the 
B-mode polarisation being perhaps the most realistic option if we are to detect it directly. The increasing accuracy 
of measurements of the CMB non-Gaussianity provides an alternative. The stress tensor of a magnetic field is non- 
linear, implying that the statistics induced on matter perturbations are intrinsically non-Gaussian, regardless of the 
nature of the underlying magnetic field. Since the standard scenario contains relatively few sources of primordial non- 
Gaussianity, it is possible that a magnetic signal is dominant. Viewed another way, predicted signals from a magnetic 
field are likely to be of a characteristic nature, and must be found in and cleaned from the CMB data before any 
conclusions on early-universe physics can be made. 



Aspects of the three-point moments have been studied in a series of papers in the last few years ("Brown &J 
Crittenden|2005t|Brown|2006}|Seshadri & Subramanian|2009l [Caprini et al.|2009[[TYivedi et al.|2010t[Cai et al.,2010'} 



Shiraishi et al. 2010b Kahniashvili & Lavrelashvili 2010 1. A bispectrum is set by three wavevectors, which we denote 



with k, p and q. Since to retain statistical isotropy these must form a closed triangle, this geometry can equivalently 
be expressed with the scalars k, r, (p, where p = rk and (p is the angle between p and q. Employing these variables the 
bispectrum geometry can be written as a foliation of planes of constant r and for each constant angle (p we then have a 
one-dimensional line through the bispectrum which in broad terms is expected to act in a similar manner to the power 
spectra. 

The magnetic bispectra studied thus far have typically been along only three such lines, all in the r = 1 plane - 
the "colinear" case where k - p - q/2 and so <p = ( [Brown & Crittenden|2005[[Brown|2 006"; Caprini et al.'2009| 
hereafter BC05, B06 and CFPR09), the "equilateral" case where k - p — q and so(p - 2;r/3 ( jSeshadri & Subramanian| 
2009t[Caprini et al.|2009[[Trivedi et al.|20T0l hereafter SS09, CFPR09 and TSSIO), and the "local" or degenerate case 



where k w p, q 
k + p. 



and so * TT (SS09, CFPR09, TSSIO). TSSIO also considered configurations where (/> = but 



The recent studies have expanded the previous results considerably. SS09 considered the equilateral and degen- 
erate lines of the bispectrum of the magnetic energy density for nearly scale-invariant magnetic fields, concluding that 
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the degenerate line provides the greatest contribution to the integral and employing an approximation to this dominant 
term to estimate the CMB signal. Likewise, CFPR09 considered the bispectrum of the energy density and considered 
the colinear, equilateral and degenerate lines. The authors generally relied on approximations that neglect angular 
terms in the integrations, or apply only on large scales. Doing so recovers the scaling behaviour of the bispectrum at 
the expense of an accurate calculation of the relative amplitudes between lines. Since the degenerate line was found 
to diverge as q^"^^ as — > this term is likely to dominate. CFPR09 also present exact solutions for the colinear case 
for both a causal field and a field relatively close to scale-invariance, which enable them to test their approximations. 
The approximations are certainly reasonable, but not ideal. In particular, since the bispectra are not positive-definite it 



is unclear whether there are strong cancellations to the degenerate line arising from other parts of the bispectrum. Cai 



et al. (2010) employed the approximations of SS09 and CFPR09 and extended the treatment to full transfer functions. 



More recently, TSSIO considered the bispectrum of the anisotropic pressure of a near scale-invariant magnetic field. 
Unlike the previous papers they evaluated the bispectrum along the degenerate hne in full, without neglecting any 
angular terms, finding it to be positive (in contrast with the approximate result, which would be negative). 

However, there are four components of the magnetic stress tensor comprising six degrees of freedom - the 
isotropic pressure/energy density, the anisotropic pressure, a vector component, and a transverse-traceless tensor com- 
ponent. Each of these are of roughly equal magnitude, as evidenced by the two-point moments ( Paoletti et al.|2009 



Brown|20T0 for example). There are many correlations that can arise between these at the three-point level, and in 
principle all of them must be considered for the CMB. The rotationally-invariant combinations are listed in BC05 
and B06, and numerical solutions for the colinear line are given, although for < -1 these took the form of noisy 
statistical realisations. |Shiraishi et al.| ( |2010c| l recently presented a formalism capable of considering the vector and 
tensor auto-correlation in full generality; our previous work considered only a fully-contracted form of the tensor auto- 



correlation and neglected the vector auto-correlation entirely. In Shiraishi et al.| ( |2010b[ l this formalism was applied 



to the 3-point vector auto-correlation. In this paper the authors employed a technique to integrate the full bispectrum 
and, therefore, evaluate an unambiguous CMB signalQ A similar paper was presented by 



Kahniashvili & Lavrelashvili 



( |2010| l. These authors focused on the symmetries of the intrinsic bispectrum and on the appearance of anisotropic and 
off-diagonal terms but in principle considered the full bispectrum. 

It is this consideration of the full bispectrum that is lacking for the other auto-correlations. This causes problems 
for causal fields and other fields far from scale-invariance in particular, where it is not necessarily to be expected 
that power is concentrated on the degenerate line. In this paper I address the full bispectrum across the range of 
{k, r, <p}. For numerical study I consider a white noise field and a field near to scale-invariance, and evaluate the three 
auto-correlations (t^Y (t^) and ^rj.^ representing the bispectra of the magnetic energy density, anisotropic pressure, 
and gravitational wave source respectively. It is also possible to find exact solutions for the degenerate line of fields 
far from scale-invariance, and I present the large-scale limits of these. (Interested readers can find the full solutions 
using the techniques presented here, or else contact the author for details.) These solutions complement two particular 
solutions presented for the colinear line by CFPR09. Due the complexity of the integration volume, also discussed in 
that work, I focus otherwise on numerical techniques. 



In section[2]l present a brief overview of the model and in section 3.1 set up the formulae necessary to find the 



bispectra. Sections [3 .2[ |3 .3| and |3 .4| consider the general case, the colinear case and the squeezed case respectively and 
then I present my results in section]?] Section ]5]provides a brief conclusion. The appendices contain some additional 
formulae and some analytical solutions for the degenerate bispectrum on large scales. 



' The authors followed this with a detailed presentation in ? while this manuscript was under review. 
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2. Tangled magnetic fields 



At the hnear level a large-scale primordial magnetic field baik, rf) decays as a^{rj). We work with the scaled field 
Bfl(k) = a^{T])ba(^, t]), which is constant over time (Subramanian & Barrow 1998a Jedamzik et al. 1998 1, and assume 
the electric fields to be negligible. Magnetic fields are damped on small scales, primarily from radiation viscosity, 
and we approximately model this effect with a time-dependent cut-off scale kdrj) above which they have no power, 
approximating a sharp damping tail. Note that there is then a time-dependence in the magnetic field Ba(x) associated 
with the damping scale which in this section we write explicitly. The damping scale freezes at photon decoupling and 



Take as a toy model Gaussian random magnetic fields with the power spectrum 



{BJk, kjT]))Bl(k', kAif))) = 'PB(k)Pat(m(kc(Jl) - k)(27Tf6(k - k') 



where 



(1) 



(2) 



projects objects onto a plane orthogonal to the wavevector k, Psik) is the magnetic power spectrum and in the interests 
of simplicity I have neglected an antisymmetric, heUcal component. H{x) is the Heaviside function, which is defined 
to vanish for < and to be unity for x>0. The power spectrum is defined as 



'PBik)^ABk"«. 



(3) 



The analyticity of VBPab requires hb > 2, corresponding to magnetic fields generated by causal processes (Durrer 
|& Caprini|2003| l. Inflation can produce fields with ub — > -3, corresponding to near scale-invariance. Some current 



constraints (generally based on the two-point moments) constrain the spectral index to < -0.12 Paoletti & Finelli 



(2010 1 at 95% confidence and ub 



-2.37^Q^3 at 1-cr. The amplitude of the power spectrum Ab can be normalised 



to observation either by the mean square field ( Finelli et al.|2008 Paoletti et al.|2009 Caprini et al.|2009 Bonvin & 
Caprini 2010 for example) or more commonly by the field smoothed on a scale A (see for example Mack et al. ( 2002| l, 



CFPR09 and BIO) but this choice is widely prevalent in the literature). In this paper we do not make contact with 
observation and leave Ab unfixed. Note, though, that the normalisation requires «b > -3 to keep the integration finite. 

The magnetic field contributes to the Euler and Einstein equations through the stress-energy tensor T{!(k). As 
the Poynting vector vanishes to first order and the magnetic energy density is equivalent to the isotropic pressure it is 
suflicient to consider the stress tensor T^(k), 



^,{^) = f\{^)5l - fl(k), where f„fc(k) = B«(k')Bi(k - k')d^k' 



(4) 



is the self-convolution of the magnetic field. The stress tensor can be separated into the isotropic pressure (scalar trace), 
anisotropic pressure (traceless scalar), vector and transverse-traceless (TT) tensor components, which we denote by t, 
Ts , tJ and r^^ respectively. For more details see B06. The magnetic field also contributes through the Lorentz force 
and this has been studied in B06 and ( [Paoletti et al.||200"9| [Paoletti & Finelli||2010| l; however, here we focus on the 
statistics of the stresses themselves. 

Magnetic fields have been separated into "ultra-violet" and "infra-red" cases by the behaviour of their stress 
power spectra; the power spectra for > -3/2 are dominated by the damping scale and referred to as "ultra-violet" 
fields, while the spectra for < -3/2 are independent of kdrj) across a wide range of scales and referred to as 
"infra-red" fields. It was noted in [Brown (20101 (hereafter BIO) that this has consequences for the "coherence" of 
the stresses. Since the damping scale is time-dependant, the statistics induced by a magnetic field on the matter 
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perturbations are in principle also time-dependant. In that paper I identified a coherence scale A^coh below which the 
stresses are "decoherent" and the standard techniques for evaluating CMB signals is suspect. The coherence scale is 
of the order of the damping scale for infra-red fields, but that for the ultra-violet fields is on much larger scales with 
^Coh ~ kc(ri)/lQQ. On larger scales the stresses then tend towards white noise. 

It is to be expected that this qualitative behaviour will hold also for the bispectra. There will be an ultra-violet 
regime in which the bispectra are dominated by the damping scale and only tend towards white noise for k < A;coh,3 
with a low A:coh,3i and an infra-red regime where the bispectra will scale as power laws. However, it is interesting 
to note that the transition between these regimes will occur not at - -3/2 but instead at ns - -1. This implies 
that magnetic fields with hb e (-3/2, -1) act decoherently with respect to the two-point moments but coherently with 
respect to the three-point moments. Extending this to higher-order moments, an a-point moment divides the ultra- 
violet and infra-red regimes at hb - -3 /a and all spectral indices in hb e (-3/2,0) can be treated coherently if one 
considers a correlation of sufficiently high order. 



3. Intrinsic Bispectra 

3.1. General Considerations 

Since a bispectrum is constructed from three wavevectors k, p and q forming a closed triangle, the bispectrum 
inhabits a three-dimensional space parameterised by three scalar quantities. I will typically take these to he k = |k|, 
r - Ipl/Aand <p = cos"' ■ p) (see Figure ill, although other parameterisations are possible (see for example Fergusson 



& Shellard ( 2007| l which employs a transformation to triangular coordinates, or SS09 and CFPR09 which employ 



{k,p,q] directly). In this paper we will work in configurations where r > 1 and <p G [0, tt]. This ensures that k < p at 
all times, while in principle q € [0, oo]. Other configurations can be found by permutations of the wavevector labels. 
This is in contrast to SS09 and CFPR09 who explicitly make the permutations clear. 

Bispectra are often evaluated in particular geometries where the three wavevectors form certain triangles, the 
most frequent of these being the equilateral case and the "local", "degenerate" or "squeezed" case in which k = -p 
and q = 0. In BC05, B06, CFPR09 and TSSIO a third situation was considered, the "colinear" case where k = p and 
so q = -2k. In terms of the coordinates {k, r, </>), these form lines of constant (p at r = 1; the colinear line is at = 0, 
the equilateral line is at (p - 2jt/3 and the degenerate line is at = n. The other configurations in the r - I plane 
describe isosceles configurations (as do certain other equivalent configurations when p = q). 

Throughout this paper I use the following nomenclature: 



"Colinear line" 


r 




1, 


/. = 


"Equilateral line" 


r 




1, 


p = 27r/3 


"Degenerate fine" 


r 




1, 


p - n 


"Colinear plane" 


r 


> 


1, 


p = 


"Squeezed plane" 


r 


> 


1, 


p — 7T 


"Bulk" 


r 


e 


[1, 


oo], e (0 



The squeezed plane is often referred to as a "local isosceles" shape, while the degenerate line is referred to as a 
"local" or simply "squeezed" shape. The colinear line has also been referred to as a "midpoint collinear" shape and 
the colinear plane as a "squeezed collinear" shape. See TSSIO for an example of such usage. 

To calculate the CMB signal we wrap the intrinsic bispectrum across transfer functions; heuristically speaking. 
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we have an integration of the form 

Bimn ~ JJf f^'^' P' 1' P' l)^Ti(k)ATm(p)^Tn(q)k^P^q^dkdpdq (5) 

with f(k, p, q, I, m, n) some function (which in general will include spherical Bessel functions and Wigner 3-j symbols). 
We therefore require in principle the full three-dimensional intrinsic bispectra, even though in SS09 and CFPR09 it 
was argued that the degenerate line dominates the integral for realistic situations. The calculations in these papers 
employed approximations neglecting parts of the integrals; as this can even cause a sign change the impact on the CMB 
signal could be significant. This was partially rectified in TSSIO where the degenerate line of the anisotropic scalar 
bispectrum was calculated exactly. However, the remaining parts of the bispectrum are still only known approximately, 



and the extent to which they modify the CMB signal is unknown. Shiraishi et al. (2010b I employed an alternative 



formalism and calculated the full integration for a vector auto-correlation. 

The intrinsic magnetic bispectra of a Gaussian field can be written (see for example BC05) as 

<T^(k)TB(q)Tc(p)) = 5(k + p + q) J r{k')'P{\V.--k:\)'P{\v + V.'\) 

xHihir^) - k'Wihiiii) - |k - \i'\)H{k, - |p + k'|) (85"ABc(©))d'k' (6) 

where [ABC] denote denote different parts of the stress-energy tensor, combined in rotationally-invariant combinations, 
and TABci®) is an angular component which is generally somewhat convoluted in form. is the set of angle cosines 
which contains 

= k ■ p, = k ■ q, 6pq = p ■ q, 
ttii^k-k', /3it=k-k-k', y*. = k ■ p -i- k', 

with equivalent expressions for Op, /3p, jp, a^, j3q and y^, and 

j6 = k' k^Hi' r = k' p'+k' /7 = k^' p'+Ti'. (8) 

In terms of the angles ^i^g and in Figure[T[ Oi^g — - cos and Opg - - cos ^j.^, with - cos 6. I focus in this 
paper on the {T(k)T(p)T(q)), {Ts(k)Ts(p)Ts(q)) and {Tjj(k)Tl\p)Tj(q)) auto-correlations. The angular- components for 
the scalar correlations are given in Appendix|Aj that for the tensor mode is extremely complex and can be found in the 
appendix of B06. A full list of the angular components of the rotationally-invariant correlations and cross-correlations, 
including correlations with the vector pieces, can also be found there[^ The formalism of Shiraishi et al. ( 2010c|a l is 
capable of dealing with the CMB bispectra of more general vector and tensor correlations. 

To tackle this equation, first write the wavevectors in units of the damping scale kc- Then if {k, r, (p) are the input 
parameters, p - rk and and 4> are the polar angles of the coordinate system, the integration mode can be written 
k' - akcirj), and 

!B = kl I a^r(a)H(l-a) I I !P (|k - a|)!P (|p -h a|) //(I - |k - a|)//(l - |p H- a|)(8nd0daida (9) 

Ja=0 •Jai,=-l J(4=() 



where k' = akr and 



|k - af - k^ + - lakau, |p + ap = p^ + + lapop. 



aj; - cos (o^ , Up - a/t cos(0) + a/^ sin(<;*) cos (^(p^ , ai; - ^ 



l-al. 



^The correlations presented in BC05 were incomplete and lacked the unwieldy tensor coiTelations. 
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This equation can be approached following a method similar to that in BIO. This is slightly complicated by the extra 
integration across 0. Note, though, that (p appears only in combination with sin(0). The cases = 0, e (0, tt) and 
= TT are then distinct from one-another and can be considered separately. This corresponds to a separation into the 
bulk, the cohnear plane and line, and the squeezed plane and degenerate line. 



3.2. 0e(O,7r) 

The integral (jojl is better cast as an integration across Up instead of (p. To do so it is useful to find an expression 
for sin 0. Writing cos (p in terms of Up gives sin = + -^a^ sin^ ^ - (Op - cos (f>)^/ai^ sin cp. Then when cp e [0, n), 
sin > and while € [n, 2n), sin < 0. 

Consider first the case where (f> e [0,n). Holding constant, dop/dcp - - -^a^ sin^ (p - (Op - cos 0)^ and so 

aP(a) I ;P(|k-a|) P (\p + a\)^-^dapda,da (10) 

Jfl=() Jat=-l Jap=a,,„ y^' 



where 



^ - al sin^ ^ - (Op - cos (11) 



and 

a prr — ak cos (f) — aii sin (f), apo — akcoscp + OksiiKf). (12) 
Now set z = Ip + a|^, transforming the integral to 

S^^f a'P(a)H(l-a) f P (\k - a\) H(l - \k - a\) f ° T ( H(l - z)^^^dzdai,da. (13) 

Here 

Zn — + + 2ap{ak cos (p — at sin zo — + + 2ap(ak cos (p + sin ^) — Zn + 4-apak sin (14) 

and 

2 2 

= ^ ■ (15) 

2p 



Considering the other case where (p e [n, 0), the derivative instead becomes dap/d(p - sin^ (p - {Op - cos (p)^ 
. Following the same process, the resulting integral is the same as in the previous case and so in general 

S=^r fl!P(fl)//(l-fl) r !P(|k-a|)//(l -|k-a|) r° !p(Vz)^(l -z)^dzdMfl. (16) 

Now, at z = Z;r and at z = zo, *P = 0. There are then poles at the limits of the integration. Between these limits, 
> 0; transforming z to z = Z;r + 2a/7/ gives ap - cos cp - atsiiKp + f and *F = f{2ak sin <p - f). is therefore 
positive. The quantity / reaches a maximum at z = Zq, where /max - '^ak sin cp. As a result, > except at z G {z„, Zo) 
where = 0. 

Finally, the integral across ak can be transformed into one across |k - a|^. Setting y = |k - a|' rapidly produces 
S=-^\ 'P(a)H(l-a) I 'P{^Jf)H{\-y) \ 'P{yrz)H{\ - z)^{W)dzdyda. (17) 

Jfl=o Jy={k-ay- Jz=z, 
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With q - -sjk'^ + + 2kp cos (p - k -y/l + + 2r cos 0, the angles that contribute to the integrations can be sum- 
marised as 

k k 

Okp - cos(p, OpQ - — ( 1 + r cos 0) , 6,, a - — {r + k cos 0) , 

q q 

k^ + - y z- - k / \ 

at- r-- , Op- , aq---[ak + rap], 

2ak zap q ^ ' 

k- aok „ k6kp - aup kO^q - aoq 

Pk = Pp = p , Pq = p , 

Vy Vy Vy 



pOkp + auk p + aup pOpq + aa, 



Jk = p , 7p = 1^, Jk 



1 



Vz Vz Vz 

- kak - a _ pUp + a _ kpOkp + a(kak - pap - a) 

P= — r= — = 1= ■ 



3.3. 



= 



When (p - Q the dependence on (p vanishes and the integral across (p becomes trivial. In this alignment aj. = 
and so 

- k^ + -lakuk, \p + = p^ + + 2apak. (18) 

These become zero at = (k^ + a^)/2ak and - -{p^ + a^)l2ap respectively. Since these occur on either side of 
Uk - Q the integration over Uk can be separated into two regions, transforming |p + a| when < and |k - a| when 
Uk > 0. Doing so yields 



S = nk 



I r aP(a)H{l-a) 
Jo 

-j^ r-k^+a- 

+ 7 'P{^|^r{\9 + ii\)H{\-yk)H{\ 

1^ J(k-a)^ 



1 np-+a- 

P{\k-a\)p(^/y;)H(l-\k-aml-yp)(Sr)dyp 



Ip + amT)Ayk 



Aa 



(19) 



where 



|k - a|2 = /t^ + fl^ + - + - y,) , |p + a|2 = + + I [k^ + a^ - yt) . 



The angle cosines in the colinear limit reduce to 



(20) 



0, 



kp 



-0 



kq 



-0 



pq 



1, Qk^ap^ -Oq, Pk = Pp = -Pq, 7k = 7p 



-71 



(21) 



and the angular terms are presented in equations ( A2 1-( A3 i. In particular, note that the tensor auto-correlation explicitly 
vanishes identically, as was argued before in BC05. 

In the colinear limit where r = 1, this simplifies further. Under the transformation yk — > yp the angular terms 
"Fttt, T^tststs and 'Ftttttt invariant even though the underlying angles change; therefore when r = 1 and = the 
bispectrum becomes 

27Tk^ r' nk^+a- I , \ 

S= ar(a)H(\-a) V {yJfjVi J2k'^ + 2a^ - y\{iT) H(l - y)H(\ + y - 2(k^ + a^))AyAa. (22) 
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3.4. -> TT 

There is a strong, intuitive difference between a squeezed bispectrum with r - I, which corresponds to the 
vanishing of q, and a squeezed bispectrum with r > 1. At r - 1, the configuration is approaching a degenerate 
state where ^ — » and at (p = n reaches k - -p, q - Q. At r > 1, however, the bispectrum is in a squeezed state 
with q - p - k. The symmetries of the bispectrum equate these squeezed configurations with an equivalent colinear 
bispectrum. Specifically, setting 

k->p, p->q, q— >k (23) 
maps !B(k, p/k,7r) onto !B{q,k/q,Q). See Figure[2] 

As a result we consider the degenerate line separately from the rest of the squeezed plane. 



3.4.1. The Degenerate Line 

Here r - I and (p - n. In SS09 and CFPR09 it was shown that for ng < -1 the bispectrum along this line is 
divergent as ^ — > with S oc q^""'^^ to leading order, so we must restrict ourselves to ^ ~1- Since o-^. = -Up one 
can write 



k^ 

k r ' 

2k Ja=0 



ar{a)H(l - a) 



[k-df- 



X2n 
( 
, =0 



(24) 



The integration across (p remains since Oq - aq{^). Here p - k and ^ = 0. The vanishing of q introduces spurious 
poles in the integrand which can be removed by Taylor-expanding around (p - n, leading to 



Okp = -1 
ait - cos 6 
k — auk 



|k-a| 

n = -Pk 

kuk - a 



|k-a| 



Okq = 0, 0pq = 0, 

Up — —at, aq — —sm6cosc 
a&k cos 4> 

Jp = -Jk, Jq = -Pq, 



(25) 



(26) 



The angular terms T are presented in equations ( A4 1-( A6 1. The degenerate bispectrum in this case is similar to 
the power spectrum integrals evaluated in (i Finelli etal.|2008; |Paoletti et al.|2009 Yamazaki et al.|2008 Brown|2()T0 



and should similarly possess exact solutions for certain values of n^. Examination of equation (24 1 shows that the 
limits on the integration are the same as for the power spectrum employed in BIO: 



For k>2, k-a>Q for all a and so 

SDeg(fc > 2) = 0. 

For k e (1,2), k + a > 1 for all a and so y^ax = l.k-a<l^a>k-l and so 

SDeg(^e(l,2))=^ r aPia) f (Piy/^ff (S^d^d^d, 

Ja=k-\ J(k-a)^ J 4>=0 



(27) 



(28) 
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• Fork e (0, 1], (k + a) < 1 

SDeg(fee(0,l]) = 



a < l-k and allows forymax = (k + a), while a> I - ^implies y„ 

l-k nik+af /-2;r 



- 1. Therefore 



2fe 



p 1 -a; /^(K+a ) r^ZTT 

Jo=0 Jy={k-a)- J(p=0 

Ja=\-k Jv=(k-af J0=O 



(29) 



Solutions exist for all integer and half-integer rig > -1/2 and are presented for large scales in appendix [B] the full 
solutions are available on request or from the author's websit^ These solutions complement the exact solutions for the 
colinear case of Sjtt presented in CFPR09, although they included a solution at «b - -2 which due to the divergence 
with q - Q cannot be found for the degenerate Une. Solutions for hb < -2 are of particular interest as these fields are 
currently the weakest constrained ( |Kahniashvili et al.|2010t[Paoletti & Finelli|2010[|Yamazaki et al.|2010| i. 



If «B < -1 this equation cannot be integrated; instead we sample finely along values of as i 
divergence and set !B{(p - n) = Oto impose numerical stability. 



7T to track the 



3.4.2. The Squeezed Plane 

Since bispectra in this configuration can be mapped onto a colinear line we expect a very similar line. (We do 
not perform this mapping explicitly as our sampling in {k, r, 0) will not cleanly transfer to a reasonable sampling in 
[k' , /, <p'). It is easier to integrate the squeezed plane directly than it is to refine the sampling to produce a suitably 
smooth result.) When r > 1 and (p - n, the dependence on drops out and so the bispectrum becomes 

Ssq = 2jikl r a^PiaWil - fl) f !P(|k - a|)f (|p + a|)//(l - |k - a|)//(l - |k - Si\){%T)dakda. (30) 

Consider the coordinate transformation |k - a|^ - y. The modulus of p + a is then 

Ip + a|2=/+fl2+^(3,,_fc2-fl2) (31) 
k 

which for fixed r > 1 reaches zero at k^ - + ryk- This only occurs at = which is at the edge of the integration 
volume; unlike the case for the colinear sheet there is no need to transform both moduli. The bispectrum is thus 



Ank^ r' r-(k+a)- 

Ssq = ^ aMa)H{\ -a) \ V{ ^Jy)V{\l^ + a|)i/(l - y)H(l - |p + a|)(8T)daidfl, (32) 

Ja=0 Jy=(k-aY 



Ankl fi . r<*+"> 

Jy={k-a)- 

with 

„2 . „2 . P 



|p + ^ + a^ + j[y-k^ - a^) . (33) 
As expected, this is similar to the colinear line. 



3.5. The Use of Statistical Realisations 

Following the techniques employed in BC05, B06 and B 10 we can imagine finding the magnetic bispectra for both 
infra-red and ultra-violet fields by analysing realisations of cosmic magnetic fields. Equation (j6| was derived assuming 



http ://f oik. uio.no/~ ibrown 
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that the magnetic fields obey an underlying Gaussian distribution; the use of realisations eases this restriction. In B06 
we analysed the colinear configuration employing a grid of side-length /dim = 192 and presented the results averaged 
from 1200 runs. These results took some time to recover and the errors rapidly grew dominant on large-scales; while 
we found good agreement between the realisations and numerical integration at A; > 0.12^^ we excised the low -A; 
region. Since it is this very region that is significant for CMB studies, this is not an ideal solution. 

Evaluating a complete bispectrum via realisations is certainly possible. However, as any general code must search 
through a grid for two complete wavemodes it is an extremely slow procedure. The scant mode-coverage on large 
scales forces one to employ the largest practical grid-size and number of realisations. Since bispectrum codes typically 
scale as A^^ ~ Z^j^^ this rapidly makes the use of realisations to generate a full 3D bispectrum unfeasible. Similar 
problems hit the attempt to employ realisations to study other one-dimensional linethroughs such as the degenerate 
or equilateral cases. In these situations the angular dependencies on force one to evaluate a subset of a full 3D 
bispectrum and the savings in time are not sufficient to render the problem readily tractable. 

We can, however, test the analysis of the colinear line via realisations. In Figures [3] and |4] the data points present 
the results from the average of 6000 realisations of a colinear magnetic bispectrum evaluated on grids of side-length 
'dim = 256. Error bars are one standard deviation. For hb - the results are in clear agreement with the numeri- 
cal integration and a significant improvement over those previously presented. For rig = -5/2 the bispectra, while 
vastly improved from those presented in BC05 and B06, exhibit a clear infra-red damping at ^ < arising from an 
unphysical infra-red cut-ofF. (This damping was discussed this damping in BIO and modelled semi-analytically with 
an IR-damped power spectrum.) For k > k^ the realisations agree extremely well with the numerical integrations. 
Unfortunately the current technical restrictions make repeating this analysis for the rest of the bispectrum unrealistic. 



4. Results 



The various bispectra, equations (17 22 19 24 and 32 1, are integrated numerically employing a Miser Monte 



Carlo routinqj 1 consider the cases riB - and = -5/2 as representative of ultra-violet and infra-red fields 
respectively. Grid dependences were tested employing a range of numbers of samples in the chains; convergence was 
found to occur at approximately n ~ 10^ for «b = and « ~ 10^ for «b = -5/2. Results were found using n = 5x10* 
for riB = and n - 5 x 10^ for iib - -5/2. The exception is the ^r^^ correlation where n - 5 x 10^ samples were 
needed for convergence across the entire plane due to the complexity of the integrand. 

The grid the bispectra are integrated on also varies with but is relatively sparse in k, particularly for - -5/2 
where the most important area is <k 1 which is expected to exhibit a scaling behaviour. This sparse sampling is 
introduced to speed the numerics; for low k the integration volume is large and integrations take quite some time, 
especially at low n^. The angle is sampled once every five degrees except near the colinear and squeezed planes 
where it is sampled once every half a degree, and around the equilateral line for ^r^^ which for hb - -5/2 contains a 
feature of particular interest. 1 consider r e [1,5]. 

Previous studies of the magnetic bispectra have typically been of scalar correlations, not least as it is only recently 
that a clear formalism for generating the CMB angular bispectra from vector or tensor perturbations has emerged 
Shiraishi et al. ( 2010c|a l. They have also been in {k, p, q} space. Employing instead {k, r, ((>} space has some benefits. 



the chief of which is the extent to which it highlights the scalings in k. In particular, the bispectra sampled along lines 
of constant {r, (/>} are expected to obey a power-law in kfoTk< kcoh- The bispectrum can then be modelled in a manner 



The Miser routine is adapted from that developed by Numerical Recipes www . nr . com. 
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similar to the power spectra in, for example, BIO as a plane at constant k and a given scaling relation. It also simplifies 
the construction of the complete bispectrum somewhat, in that we can integrate cleanly along lines of constant {r, 0) 
and assemble these lines into a 3D bispectrum. The other main benefit of employing this coordinate system is that 
it results in a fine sampling around the r = 1 plane; from the results of SS09 and CFPR09 and the decaying of the 
bispectrum as r — > oo, the degenerate line should dominate and other configurations in and close to the r - I plane 
will be more significant than those further away. 

To get the full bispectrum in {k,p,q} space one can transform the results. Since they obey the triangle relations, 
bispectra in {k, p, q} space are wedge-shaped. Transforming to this space yields a cone filling one third of the wedge. 
The other two thirds of the bispectrum in {k, p, q] space can be found from permuting across k, p and q. 

There are, however, drawbacks to the coordinates {k, r, <^}, the chief being a lack of immediate clarity. In particular, 
the q^""^^ scaling found near to the degenerate line for rig < -1 in CFPR09 is obscured. Further, reconstructing the 
bispectrum in {k, p, q} space will be slightly complicated when considering cross-correlations such as (T(k)T(p)Ts (q)) 
since the anisotropy along the q vector must also be taken into account when performing the permutations. However, 
I feel that the convenience with which the bispectra can be found and, more importantly, the clean scalings in k which 
should be recovered outweigh these disadvantages. 

The bispectra are assembled with the following procedure: 

• For a set {r, (p) sample a bispectrum along a line in k. 

• Repeat this for a set r and samples in € [0, n]. 

• Assemble these linethroughs into a plane at constant r, interpolating between missing points when necessary. 

• Repeat this for samples in r € [1,5]. 

• Assemble the planes of constant r into a 3D bispectrum. 



There are a few general statements can be made about the form of a bispectrum before calculating the results. 
Along the colinear line the bispectrum will have power until k = 1, which can be seen from evaluating the limits 



in equation (22 1. Conversely, along the degenerate line, the bispectrum will have power until A; = 2, as seen from 
equation (27i. We expect the power in the r - 1 plane to smoothly connect the two. As a triangle on the squeezed 
plane can be rotated onto the colinear plane these parts of the bispectrum should be expected to resemble one-another, 
except when r = 1 and the rotation breaks down. 

In the same way, the limit r — > oo corresponds to the degenerate line on small scales since large r implies that 
k p and k q. The degenerate line itself will have a coherence length of A:coh ^ 1 and on smaller scales is 
decoherent. In the limit r — > oo, then, the bispectrum should be expected to be decoherent across almost all scales. 
The coherence scale therefore obeys ^coh('" ^ oo) <s; 1 . The conclusion is that the coherence length is a function of r 
and decays for high r, although it is not necessarily a function of 0. 

Furthermore, since for near scale-invariant fields one expects S ~ k^"B+^ across at least most of the bispectrum, 
the power on small scales is dramatically less than that on the largest scales. Another corollary of the equivalence of 
the high-r regime and the degenerate line is then that the power in the bispectrum should decay as r — » oo. It will 
also damp entirely when p > 2 which occurs at ever-smaller k for increasing r. As r increases, then, the bispectra for 
infra-red fields should lie within a narrowing wedge, decaying in amplitude at the same time. At higher r the signal 
will be decoherent across almost the entire range. For such fields the decoherence at high r is unlikely to be significant 
on the CMB. 
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For ultra-violet fields this is not the case since the large-scale behaviour is expected to be white noise, and the 
small and large scale signals are of approximately the same order of magnitude. While the coherence length will still 
evolve to larger scales as r increases and the bispectrum damp away when p > 2, the amplitude will remain of the 
same order-of-magnitude and the decoherence at high r will be significant. This implies that our coordinates are not 
well-suited for the study of ultra-violet magnetic fields on the CMB. On the other hand, the caveats mentioned earlier 
and in B06 and BIO concerning the validity of wrapping decoherent statistics onto the CMB suggest we should take 
care with such cases regardless of the coordinate system employed. 

The general structure of the bispectra should then be clear. There is power in the r - I plane ranging from k - kc 
at the colinear line tok - 2kc at the degenerate. For increasing r the power dies at increasingly small scales until in the 
Umit r ^ oo it tends towards the small-scale (large-fe) degenerate solution. The coherence length A;coh likewise decays 
from a maximum in the r = 1 plane until it vanishes at r — > oo. The squeezed and colinear planes exhibit a similar 
behaviour to one-another although the quantitative nature is somewhat different. For infra-red fields the amplitude at 
high r will be negligible, but for ultra-violet fields it will be comparable to that on the r - I plane. 

In analogy with the two-point moments of the stresses, on large scales a complete bispectrum for a typical 
magnetic power spectrum (such as the power-law considered here) will take the form 

m,r,<f>)^S.(r,4,)lJ-^ (34) 

for some pivot I. The values a - for ng > -1 and a = 3??^ -H 3 for < -1 are expected, although this 

must be verified. This reduces the CMB integration to an integral across transfer functions and only a 2D plane rather 
than a 3D plane. In some situations it is possible that analytical approximations might be found given the power-law 
behaviour of the bispectrum. In any event, the full integration is then reduced to finding S*(r, 0), the full knowledge 
of which is the factor missing in previous calculations. (Indeed, when TSSIO evaluated the degenerate line in full for 
an infra-red field they found that the full result has an opposite sign to the approximation neglecting S*(l, tt).) It is far 
faster to integrate a plane with <p e [0, tt] and r e [1, r^ax] for constant k than it is to integrate a full bispectrum. Such 
an increase in speed is vital if we are to perform accurate parameter estimation. 

In principle the plane S*(r, <p) should be tilted such that fe*(r) < fecoh('"); or chosen to be so small that A;coh('") 
only grows through A:* when the power in the bispectrum has become negligible. In practice, for infra-red fields the 
degenerate line is expected to dominate the bispectrum so significantly (SS09 and CFPR09, for example) that it is 
unlikely a significant error will be introduced by selecting a constant A:* » (3(10"'* - 10"^). 

In the results that follow, the errors presented are 1-cr. These errors are dependant on the convergence of the 
integral as well as the number of samples taken in the average, owing to the Monte-Carlo nature of the integration. 



4.1. The Flat Spectrum, = 

Consider first the stresses of a white-noise magnetic field. This case is straightforward to integrate as it contains 
few of the poles that litter the integrations for infra-red fields. Analytical solutions are, however, extremely difficult to 
find and it is quite possible that they don't exist, with the exception of solutions along the degenerate line (see Appendix 
[B) although some solutions can be found in the small-scale tail. I integrate the equations numerically, setting a grid- 
spacing in (p of five degrees and employing 5 x 10^ samples in the Monte-Carlo chains. It is to be expected that the 
bispectra tend towards white noise (a - 0) for k ^ I with a coherence length on the order of kcoh ~ IQ-^ - 10-^ in the 
r = 1 plane. The magnitudes of the bispectra should be of the order of (1 - 10)Ag, although it is impossible to predict 
what configuration, if any, dominates. 
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4.7.7. Srr. 

The auto-correlation of the isotropic pressure is straightforward. It is positive definite, showing a clean evolution 
between = and (p - n and between r - \ and r - 5. Figure [3] shows the colinear, equilateral and degenerate 
configurations. The solution for the colinear line agrees well with the results from realisations which are, unfortunately, 
extremely noisy at low k, albeit a great improvement over those previously presented. 

Figure |5]presents slices through the bispectrum at r e { 1 , 2, 3, 4). In the plot for r - \ the transition between the 
colinear and degenerate behaviours is clear, with power decaying on ever-larger scales with increasing (p showing only 
a minor evolution for increasing r. The simplicity of Tjjj implies that integration is effectively constrained purely by 
the Heaviside functions in the power spectrum and not by cancellations in the integrands. The volume this bispectrum 
fills can then be taken as approximately marking the extremal bounds of any of the bispectra - that is, power should 
not be expected outside of the bounds of (t^^ Furthermore, the amplitude of the white-noise signal is almost constant 
with both r and (p. This is in agreement with and extends the results of CFPR09 who found that, for ub > -1, the white 
noise amplitude is constant across the colinear, equilateral and squeezed configurations. 



Figure 11 gives the (r, 0) plane at a pivot scale A;* - 7.5 x 10"^, which shows this strikingly. Modelling the 
bispectrum with a power-law around this plane, and for presentation choosing to evaluate a at the equilateral line, the 
integrations here give 

Q- = (0.16 ± 3.23) xlO""^ (35) 

This is the result of sampling the average of six points in the white-noise regime. Considering the full plane a(r, <p) at 
k - ki, and low r gives results which are similar The coherence length evaluated at the equilateral line A;coh evolves 
slowly towards larger scales from kcohir = 1) =s 5 x lO^^* as r increases. A regime with A;coh('") ^ ^* is reached at 
larger r. The recovered a would grow inaccurate unless = A;*(r), reducing at large r. The constant back plane and 
the white-noise behaviour of the scaling make this bispectrum particularly simple on large scales. 

Approximating S*(r, (p) - 8.35 and a -Q, then, with a pivot fe*(r) that always remains within ^coh('") produces a 
good model of the bispectrum as it is applicable to CMB scales. However, the rapid decay of fccoh with r, combined 
with the relative strength of power on small scales compared to large scales, implies that one should be wary of 
wrapping these fields directly onto the CMB. 



Finally the full bispectrum in k,r,4> space is presented in Figure 13 This plot shows isosurfaces of constant 
amplitude through the bispectrum. This figure highlights the simple sheet-like structure of this correlation, the wedge- 
shape of the bounds of the correlations and the relatively slow decay with increasing r. 



4.1.2. S-r^rsTs 

The auto-correlation of the anisotropic is more interesting in structure than that of the isotropic pressure. In 
amplitude it is similar, and on white-noise scales is Aj^ » 0.91 At, but it has a much more complicated pattern. Figure 
[3] shows that the signal is negative along the entire colinear line, while it is positive for much of the range of k along 
the equilateral and degenerate lines. Even here, though, it is negative on small scales. The ^t^^ signal therefore passes 
between regions of positive and negative power dependant on both <p and k. It should be expected that the same is 
true for varying r and this is confirmed in Figure |6] which shows the negative region beside the colinear line, which is 
joined by a corresponding region along the squeezed plane. As expected, away from the degenerate line the squeezed 
plane has the same general behaviour as the colinear plane. The planes in Figure |6] slowly squeeze a positive region 
which is centred approximately around = n/2. The coherence length at the equilateral line is A;coh ~ 10"^, again 
evolving only very slowly with increasing r. 
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From the general arguments above, in the limit of high r we would expect the bispectrum to tend to the large-scale 
degenerate solution. In Figure [3] we can see that this is negative, albeit extremely small. The central, positive, region 
is then expected to entirely vanish at the extreme reach of the bispectrum. As with the isotropic pressure, the signal 
along the colinear line is well-modelled with statistical realisations, but is extremely noisy on large scales. 



Figure 1 1 shows S*(r, 0) at fe* - 7.5 x 10"^ and demonstrates these features clearly. The region of positivity is 
asymmetric around the equilateral line at = 2nl3 with a tail at r = 1 extending to reach the degenerate line atcp-n. 
The bispectrum in general exhibits a characteristic "triple drainpipe" pattern. Recovering the scaling relation around 
the equilateral line gives 

a = (-1.73±4.97)xl0"^ (36) 
The bispectrum can then be modelled by using S*(r, 0) and a = 0; unfortunately this plane is not immediately 



amenable to analytic approximation. The full bispectrum in Figure 13 shows the triple-drainpipe extremely clearly, 
along with its slow decay at high r. 



4.1.3. ^Bfj-Tj-Tf 

As with the isotropic pressure, the bispectrum of the tensor anisotropic stress is positive-definite, with a peak 
magnitude of A-,-j. ^ 0.48A^. However, unlike the isotropic pressure, this signal has a non-trivial structure. We 
commented in BC05 that one can quickly prove that the colinear line of this signal is identically zero by the symmetries 
of the projection tensor, a conclusion which the full numerical integration confirms in Figure [3] which shows the 
vanishing colinear signal, but the equilateral and degenerate lines being of the same order-of-magnitude as the scalar 
auto-correlations. 

Indeed, examining the plane at r = 1 in Figure |7] reveals that the result is negligible (although, in principle, 
non-vanishing) up until (p ^ n/S. For higher r the symmetries between the squeezed and colinear planes imply that 
the squeezed plane should also vanish, and this is observed. By r = 4 the signal vanishes for < n/4 and (p > An/S. 
The maximum wavenumber at which a signal is non-vanishing is also strongly dependent on <p, reaching its peak at 
<p ~ 0.5771 for low r, diminishing slowly with increasing r. This selects a configuration of particular interest - an 
isosceles lying between an equilateral and a right-angled triangle - for this signal which could not easily be predicted 
from examining the angular integrand. 

The coherence length at the equilateral line is kcoh ~ 5 x 10"^ for r = 1 and as with the scalar modes evolves only 
slowly. 



The full bispectrum shown in Figure 13 reveals that, as r increases, the region with a significant signal continues 
to narrow. It also highlights the characteristic shape of this signal, similar to a tree trunk, with a bulge on the base 
stretching out to the degenerate line. Sj^tttt^ plotted in Figure 11 shows that the power on large scales in this 



signal, in contrast to the two scalar auto-correlations, is focused in the vicinity of the r - 1 plane and near-equilateral 
configurations. The high-r solution is, however, similar in order-of-magnitude to the degenerate hne, as should be 
expected. The scahng recovered around the equilateral line is 

ff^ (-1.52 ±5.17) X 10"^ (37) 



The bispectrum can then be modelled with the plane S*(r, (/>) from Figure 1 1 and a - 0. 
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4.2. The rated Spectrum, = -5/2 

Consider now the stresses of the ultra-violet field with hb - -5/2. These integrals are much more difficult to 
integrate since the integration volumes are steeply tilted towards poles at the edges of the integration. While such inte- 
grals can be controlled numerically, much longer chains are required if convergence is to be reached: 5x10^ samples 
for the scalar modes, and 5x10^ for the tensor auto-correlation, which possesses a significantly more complicated 
integration surface. Analytical solutions are difficult to find and may well not exist. As with the white-noise field, is 
sampled once every five degrees, except around the equilateral line for ^t^ ^ which contains an interesting feature and 
is sampled at every half a degree. 

The bispectra are expected to scale along lines of constant {0, r) as S oc A:^*""^''. This scaling has been observed 
for colinear (BC05, B06, CFPR09) and equilateral (CFPR09) bispectra and is naively expected to hold throughout the 
rest of the bulk. Towards the degenerate line, SS09 and CFPR09 found divergences S oc q^>'n+^ and concluded that 
this line dominates. In the {k, r, 0) coordinates this scaling will be obscured. Specifically, the dominant term found in 
CFPR09 was ^ 

S - 5(k + p + q)-^ f , , — -rXq'"'^' + ...); (38) 



U4n^ \(nB + 3)(2nB + 3) 
after transformation this becomes 

« » * - p - ( „,..3),;,..3) (■ - ^f""'" - 4 

The scaling in k would then be expected to be the same as for the other configurations; however, a large drop in 
amplitude is expected with only a slight change in r and (p. For the field with ub - -5/2 the scaling should be 
a = 3{nB -H 1) = -9/2. 

The stress power spectra on small scales for infra-red fields, k > 1, are of the same order of magnitude as those 
for ultra-violet fields (see, for example, BIO). This behaviour is expected to hold for the bispectra, too. Given the 
scaling on scales k < kcou, very large signals in the regime applicable to CMB studies are then expected. 

In principle, the bispectra should lie within the same volume occupied by SttAub = 0). The strong scalings, and 
the expected domination of the degenerate line, will most likely obscure this feature. 



4.2.7. Srrr 

The ID bispectra along the colinear, equilateral and a near-degenerate line are shown in Figure |4] compared, 
in the colinear case, with the results from realisations. The agreement with the realisations is reasonable but the 
consequences of the infra-red cutoff" are evident with the evaluation losing power on increasingly large scales. The 
signal is positive-definite, which agrees with the signs found in CFPR09 and SS09. Interestingly, the equilateral line 
is of the same order-of-magnitude as the colinear - specifically, Seq/Scoi ~ 1-02. This is in conflict with the results 
of CFPR09 where the authors found that Seq/Scoi ~ 0.1. The discrepancy is almost certainly due to our inclusion of 
angular terms necessarily neglected in the approximations of CFPR09. In contrast, the bispectrum evaluated along a 
line at = 357r/36, corresponding to = 175°, is three orders of magnitude greater Lines nearer degeneracy are, 
naturally, larger again. As expected, the degenerate line will dominate the integration. 

Figure[8]shows the bispectrum sliced at r e {1, 2, 3, 4). These plots show the evolution of the bispectra from <p = 
Xo (p - 71 clearly and the slowly-diminishing power as r increases. The most dramatic feature is the area near to the 
degenerate line where the power diverges. By r = 2 the squeezed plane has approximately equal power to the bulk. 
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The lower curves in Figure]?] damped at small k, highlight this; these show the colinear, "equilateral" and squeezed 
configurations for r = 5. While the colinear and "equilateral" signals are only approximately an order-of-magnitude 
lower than that at r = 1, the squeezed signal is suppressed by five orders of magnitude. 

The coherence length fccoh evolves from A;coh ~ 1 at r = 1 to kcoh ~ l/10atr = 5 and therefore generally remains 
on much smaller scales than the CMB. 



Figure 12 shows S*(r, 0) at a pivot scale of 4.4 x 10"^, where the pivot is taken to be larger pivot than for tig = Q 
since A;coh is reached at a smaller k. The scale has been truncated to show the increased structure in the pivot plane 
compared to the ng - case. As with those ^r^^ and (r^'j cases, this is unfortunately not immediately amenable to 
analytic approximation. Modelling the bispectrum with the power law around this plane we recover 

a = -4.498 ± 0.006 (40) 

which is extremely close to the predicted a = -9/2. Using the S*(r, <p) plane with this value of a can then recover the 
section of the bispectrum relevant to the CMB. 

The extent to which the power piles onto the degenerate line renders diagrams of the full bispectrum at - -5/2 
somewhat unhelpful. The most useful information can be extracted from the lane presented in Figure ]T2] and the 
recovered value of a. It can also be noted that the increase of the error in a can be used to diagnose when fccoh('") has 
grown through fe*(r); however, given the domination of the degenerate line it is not expected that much error will be 
introduced at high r from employing a constant fc*. 



4 2 2 !R 

As with the previous bispectra (both for ng - Q and ng - -5/2), the r - \ plane can contain features that only 
exist right on this plane due to the high symmetries that can occur when k - p, the clearest examples being that the 
degenerate line for = is positive for ^r^^ and non-vanishing for (t^), in contrast with the behaviour in the rest 
of the squeezed plane. The ^t^^ bispectrum at rig - -5/2 illustrates this even more dramatically, although it also 
contains other unique features. 

Figure |4j shows the colinear, equilateral and near-degenerate lines. As with ^t^^ the agreement with realisations 
for the colinear line is good, with the reaUsations exhibiting an infra-red damping. The two scalar signals are indis- 
tinguishable from one-another along this line. Along the degenerate line, the signal has become negative, and almost 
equal in magnitude to {r^')- This is in agreement with TSSIO, recalling that their Ylg - -ts, who found that their 
degenerate line was positive in contrast to the approximate solutions. 

The equilateral line is more interesting. The signal along this line is significantly suppressed across much of the 
range of k before flattening out as A; — » 1. For ^ > 1 it is of the same order of magnitude as the other two correlations. 
The bispectrum of the anisotropic pressure then exhibits some interesting behaviour in the r = 1 plane, passing from 
positive to a negative region near the equilateral line. It remains negative near the degenerate line. 

However, the linethroughs at r = 5 in Figure]?] and the planes in Figure ]9] demonstrate that the signal does not 
remain negative between the equilateral and degenerate lines. Instead, there is a negative trough, or "river" running 
between regions of positivity before a strong negative divergence as <p ^ n. Conversely, the small-scale signal is 
negative for almost all (p. The manner in which this bispectrum changes signs on large scales could potentially lead 
to very characteristic signatures on the CMB. Furthermore, the fact that the degenerate line possesses a different sign 
to the scalar signal implies that the total magnetic signal - the sum of all possible CMB bispectra, covering auto- and 
cross-correlations of all the components of the stress tensor - will exhibit cancellations. 
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The planes at r e {2, 3,4) demonstrate further oscillations in the signal; a negative region at A: — > A:^ grows in area 
and then declines again with increasing r. However, these interesting "caves" in the bispectrum are unlikely to leave 
significant signs on the CMB as the transfer functions and the amplitude of the bispectra close to the damping scale 
are negUgible compared to the largest scales. 

In TSSIO the authors also evaluated the result for a "local isosceles" shape, which they found to vanish. No such 
signal is seen here, although naturally some regions of the bispectrum lying between regions of positive and negative 
value are zero. It is likely that the inclusion of terms they neglected to enable an analytical integration account for the 
discrepancy. 



Evaluating the pivot plane at ki, - 4.4 x 10 , presented in Figure 12 and the scaling at the equilateral line yields 



a = -4.52 ±0.12. (41) 

This value has been recovered sampling only four points - the levelling of this curve seen in Figure |4]contaminates the 
calculation otherwise. The downside is the increase in the value of the error, which may be mere statistical fluctuation 
or may be a hint that the integrals should be run at a slightly higher resolution although convergence is still reasonable. 
Repeating the evaluation along the colinear line (which is very well sampled since there is one less layer of integration) 
yields 

a = -4.4999 ± 0.0003 (42) 

which dramatically confirms the scaling to very high precision. S(r, 0) shows the pivot plane, plotted on a truncated 
scale to highlight structure, is non-trivial in nature. There are negative regions at r = 1 with <p - In/ 3 and (p ^ tt; the 
inset shows a close-up on the region around the degenerate line. The coherence length behaves as for the S^^- 



4.2.3. 



Finally, consider the (tj^ signal. The colinear, equilateral and degenerate lines in Figure|4jshow that the colinear 
signal is identically zero while the equilateral and degenerate lines are of the same order of magnitude as those for 
{t^Y with the tensor signal slightly dominating the equilateral line and the scalar signal dominating the degenerate 
signal. As with the other auto-correlations the degenerate signal is strongly dominant. The signal is positive. 



The slices at constant r in Figure 10 reveal a more complicated structure. The colinear plane is identically zero 



and, by rotation, this implies that away from the degenerate line the squeezed plane is also zero. A negative region lies 

4.4 X 10"^ 



next to each of the planes. Figure 



12 



shows the evolution of these regions as r increases at a pivot of ki. 
Otherwise the signal is positive-definite. The scahng recovered along the equilateral fine around this plane is found to 
be 

a = -4.500 + 0.003. (43) 



The pivot plane S*(r, <^) from Figure 12 can then be employed with a 
bispectra, non-trivial. The coherence length again behaves as for Sji . 



-9/2. This plane is, as with most of the other 



5. Conclusions 

This paper presents numerical integrations of the auto-correlations Sjtt, ^tststs ™d Stjtjtt which are required 
for an evaluation of the CMB angular bispectrum imprinted by magnetically-induced scalar or tensor modes. On 
large scales, quantified as lying below a coherence length kcoh - kcoh(nB,f',<p), previous authors have found that the 
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bispectxa from ultra-violet fields with ^ -1 act as white noise, while those from infra-red fields with < -1 scale 
with k^""'^^, except along the degenerate line where they scale as q^""^^. This behaviour is confirmed, ultra-violet fields 
having a coherence length kcoh ~ kdrf) /IQQ and infra-red fields having a coherence scale k ~ kdrj), where the values 
are taken at r = 1 . 

In the region applicable to CMB studies the bispectra are well-modelled by a 2-dimensional plane S*(r, <p) sam- 
pled at a pivot of A: = fc*. The ultra-violet bispectra typically have power up to a large r, with A;coh decaying with 
increasing r, so if one were to put an > ~1 bispectrum onto the CMB in principle one should set a pivot on very 
large scales and allow it to run with r. In practice this would complicate the CMB integration significantly and, while it 
is not clear that the standard formalism directly applies to such fields, care should be taken. For hb < -I the evolution 
of ^coh is less important since this remains approximately on the order of A;coh ~ 0(1) except in regions when the 
bispectrum is negligible. 

These 2-dimensional planes are non-trivial but (relatively) quick to evaluate since evaluation is only needed at a 
constant A:*. The sampling can be made arbitrarily fine as r — > 1 and (p ^ n to capture the most important region, 
and can be relatively sparse outside of these regions. To constrain magnetic parameters from the CMB parameter 
estimation is required, which involves many multiple evaluations in a Monte-Carlo chain, so this increase in speed 
is necessary. There are two possibilities: evaluating S*(r, 0) as requested for each value of ng, or pre-evaluating a 
number of planes with discrete choices of «b and interpolating between them as required. The first case would be 
the ideal but might unfortunately remain prohibitively slow. The latter option rests on the assumption that the plane 
S*(r, 0) evolves smoothly with varying ng. In this case one would imagine sampling choices of hb e (-1, -3) and 
stacking these planes into a cube in {r, <p,n shspace. This issue is under current study. 

The CMB signal from the scalar modes is being increasingly studied (see for example SS09, CFPR09, TSSIO 
and Cai et al. ( 20 1 1) and that from the vector auto-correlation has recently also been evaluated ( Shiraishi et al.|2010b 



|Kahniashvili & Lavrelashvili|2010j ). An immediate consequence of this study is that we can evaluate the CMB angular 
bispectrum arising from magnetic gravitational waves and we will present the results in a forthcoming work. 

Necessary extensions to enable a full parameter estimation include evaluating the scalar cross-correlations > 
Sttsts^ ^tststs7 the scalar/vector and scalar/tensor cross-correlations SrT,,Tv^ ^tstvtvi ^ttttt^ ^ttttt ^i^d 'he vec- 
tor/tensor cross-correlation SryTrTv^ the colinear lines of all of which were considered in BC05 and B06. In the 
first instance we would want to select a constant spectral index hb ~ -5/2 and evaluate the bispectra across a section 
of the coherent range to confirm the nature of the scalings and the behaviour of permutations of the wavenumbers. For 
more general it should suffice to evaluate instead the plane S*(r, 0) at ki, * 10 We may however find for CMB 
integration that, given the form of the integration ( Ferreira et al.|1998 Wang & Kamionkowski|2000 Shiraishi et al. 



2010c I swapping back to {k, p,q]-space is more convenient. This would depend on the balance between the speed- 
increase from considering only 2D planes and the complications of remapping from the {k, r, 0) coordinate system, 
incorporating as well the permutations for cross-correlations, and is currently under study. 

Since BC05 and B06 demonstrate that each of these cross-correlations is non-vanishing and of the same order- 
of-magnitude along the colinear line each one will, in principle, leave traces on the CMB. Furthermore, the fact 
that the degenerate line possesses a different sign to the scalar signal strongly suggests that the total magnetic signal 

- the sum of all possible CMB bispectra, covering auto- and cross-correlations of all the components of the stress 
tensor - will exhibit cancellations. In particular, since magnetic bispectra are not positive-definite, contributions from 
unexpected features such as the "river" in ^t^ ^, could introduce cancellations in the integrations particularly for indices 
farther from scale-invariance where the degenerate line isn't so dominant. Perhaps more importantly, the complete 
set (Stt-t^ 2^ttts V ■ ■ ,STyTTTv^- ' ' ) contrlbutcs to the magnetised CMB angular bispectrum and the nature of these terms 

- particularly their sign near to the degenerate line - is unknown. Whether these efi'ects have a significant impact 
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on the CMB bounds is uknown at present, but conceivably they could significantly tighten (or weaken) them. Until 
the complete set of intrinsic bispectra are known, and then wrapped onto the CMB with the corresponding transfer 
functions, CMB constraints on magnetic fields from bispectra should be treated with some caution. 



Many instructive and useful conversations have helped me while I was undertaking this study but I owe a partic- 
ular debt to Kishore Ananda, Chris Byrnes, Frode Hansen and particularly Richard W. Brown and David Parkinson, 
who assisted with some of the codes employed in the project, and Robert Crittenden who introduced me to this topic. 



A. Appendix A: Angular Terms 

The angular terms in the bispectra for the scalar auto-correlations are given by 



6 



Wrrr ^ P + f+f-m and Tr.rsrs = ^ '^"^^srs (^1) 



n=0 



with 



- 87"° = 9 



TsTsTs 



-^"^TsTsTs = S^JiiPkn + PpTp + Pqjq + ^Pr) + jiam + apy,, + aqyq) + ~P(,akPk + OpPp + aqPq) 

+30kpiakap +PkPp + ykYp) + 30kq(akag + PkPq + ykjq) + 36pg(apaq +PpPq + ypyq) + 96kp9kqdpqj, 

-STr.TsTs = -3(yjlakPk +PJiapyp +PyPqyq + 3QJ9kpPkyp + yOkqakJq +p0pqapPg) 
+3{akPk(apPp + aqPq) + apyp(akyk + a^r?) + PqJqiPkyk +PpJp)) 

+9(Okp0kq7prq + OkpOpqPkPq + G kq^ pqU kO p)^, 

-^'^TsTsTs = ^{jlakPkapyp + yakPkPqyq + PapypP^y^ + 3{ekpPkypPqyq + OkqOkapypyq + OpqOkPkapPqij, 
-8r,]rsrs = -27akPkapypPqy,. 
The angular term for the tensor auto-correlation is extremely unwieldy and may be found in B06. 
In the colinear limit Ttststs ^"d Ttjtttt reduce to 

„^ q/ 1 2^r,2^ 2 ^f' +Ji^ -PW ^akykM + akPkW+PkJkPy 
^TrsTsTs = 9 - 1 + Q-j +p^ +y^ + + 



9 3 
-"IPl - "hi - Phi - akPkllP - "kPhkJ - alPkJkJi + MPhi) (A2) 

and 

(A3) 



where in the second equation we have used the definitions of the angles given in equation (21 
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In the degenerate limit when = ;r the angular terms become 



= 1 (A4) 



= 1 - 3 (a^ + ^2 ^ + - - yS^) + ;8 +3 (at (n - A) - 2a,y8,)ye (A5) 

and 

9>rTrrr = -3 + 3Qr^ + 4jS^ + 6^^ + + 3a^ + 5j8^ + 3jS^ - 3 (a;t06i - y^) + la^Pq)'^ + 3a;t 08* - jk) a^fig 

+4 {Pl+yl + 3y8m -yS^) (27,' - - -f] - yl [ic^^ +13] +f) (A6) 

-Pkjk [a] + 2(5] + 2fY pl (a] + pl+fy akih - TkWlP + + Jkf a^p^p 
+2agPlp + al {2PkJk [Pkjk - pf) +P^)-ak(Pk- Jk) a^p] - Pkjkalpl 



B. Appendix B: Ultra- Violet Degenerate Bispectra on Large Scales 

Analytic solutions can be found for the degenerate bispectra across the entire range of A: for ng > -1 and are 
available from the author on request or on his website. For simplicity only the Laurent expansions around = are 
presented here, representing the degenerate bispectra on the largest scales. Large-scale solutions in this regime can 
also be found for the coUnear bispectra and the general case by reducing the angular terms to the large-scale limits, 
which depend only on aj^; in these cases the expansions here are accurate only to hnear order. 

In the infra-red regime when hb < -1 the solutions include a term dependant on the logarithm of an infra-red 
cut-off scale. For example, on large scales when hb = -3/2 we find 

1116 32 8 16 \ 16 ,\ , 

Srrri<P = = 1 1 ^ + Y l"^^) " 3 + y ln(^/Un) 1 A"'^' " y " " nk] (B 1) 

where kmm is the infra-red cut-off. Since there is no immediate model motivating such an infra-red cut-off (barring a 
realisation on a finite grid) these solutions are not presented here, although they are straightforward to recover. 

The results presented in this section complement those in Appendix E.2 of B06 and the more general but colinear 
results of CFPR09, extending them to the degenerate bispectrum and the anisotropic and tensor auto-correlations. In 
both cases we extend the solutions to a wider range of hb- Comparison with the full solutions for > -1 demonstrates 
that these expansions (to third order in k) are vahd for k < kc{ri)l2Q. 



When Hb = 3, 



-^Trr — ^3 .iA. -t- 5 A. 3 /c j'^K^, x>tjtsTs - 105 8 75 240 J'^'^c' 

<R - 68_ _ 191, , SMp _ 107^3^_p 

■Otttttt - \3i5 24" ^ 525" 45 " )"'^c- 



(B2) 



WhennB = 5/2, 



« - ("16 _2ik-l- 1221,2 _ 59p\ p a, - ('272 _ 1- _ 290p , 509t3\ p 

« _ C 544 19/ 852.2 4.781 , 3\ .3 

^tttttt 1,2,205 24 ^ 595 2,880 /f 



(B3) 
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When Hb = 2, 



■Ottt — y^ 21'' 3" ^''''c' -^T-jTsTs - (^315 105" ^ 24" ^''''c' 

<H 272 _ 191, , 124.2 _ 77 

■^tttttt — \g45 24 105 72 )"'^c- 



(B4) 



WhennB = 3/2, 

« - fie _2fc+ ^2fc2 _ lit3\ p g - C272 _ l-^_ _998_t2 , 55p\ p 

« - /"J^^, _ 19r. , MZ6I.2 _ 349p\ p 

■^TTTTTT - (^1,575 24'*'^ 1,155^ 576^ ^'^'e- 



Whenng = 1, 

« - C 136 _ lit . 24^2 _ 9It3^ ^3 

■^ttTtTt - \3i5 24^^35'' 360" ^''"c 



Whennz, = 1/2, 

« - fie _2it- ijt2 + 13t3\ 1^3 g - f 272 _ 1,_ 706.2 , 137. 3^ ^3 

■O-m-yg 24" ^-^"c ■'-'tststs - \m g'^ 525" ^ 480"^''"c' /r>7\ 

« - C544 _ 191. . 236p_ J73_p\ p 

■"ttTtTt — (^945 24" ^ 525" 2,880" ^''"c- 



When Mb = 0, 



S^rr = {l-2k- Ik^ + ^ l)k^)nkl S,,, = (1 - - fffc^ + (i + i;r»^3, 

■Dtttttt - 24"^ 15" ^ Us 64'^ j" j'^'^c- 



(B8) 



Whenns = -1/2, 



« - ("IS _2ik- 128i;3/2 , 17,2 , i3,3\ ,3 g _ ('272 1, 79,888 , 3/2 , 82 ,2 . 19/,3\ / 3 

« _('544 19., 608 p/2 12.2 13 

■"tttttt — \3i5 24" ^ 29,835" 35" 576" ^""c' 



(B9) 



Finally, when hb = - I, 



Srrr = (4 - 81n(A;) - 2/t + f fc^ + Ik^) nk], Sr.r.rs = (-fll " W " + li*^ + Ig'^^)^'^?. 

« _ ('39,472 272 . 19, 44 ,2 7 , 3\ / 3 

■DTT-Tj-Tr - 11025 105 24" 105" 36" ^'^"c- 



(BIO) 
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Fig. 1.— 



The geometry of a statistically homogeneous 3-point moment and the integration wavenumber k'. 



P 



a) 
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b) 



c) 



d) 



Fig. 2. — Bispectra (f> - and (p - n. a) The colinear configuration, b) The degenerate configuration, c) (p - Q. d) 
(p - ji. Note that case d) can be transferred into case c) if we map q — > k, k — > p, p — > q. 




Fig. 3. — The colinear (0 = 0, left), equilateral {(p - 27r/3, middle) and degenerate ((p - n) limits of the (^t^^ (black), 
^T^^ (blue) and ^r^^ (green) bispectra, for iib = 0. 



Fig. 4. — The colinear {(p - 0, left), equilateral {(p - Injlt, middle) and near-degenerate (0 = QSlln) limits of the 
{t^"^ (black), ^T^^ (blue) and ^r^^ (green) bispectra, for hb - -5/2. Lines decaying for A: > 1 are at r = 1 while those 
decaying at smaller ^ are at r = 5. Dashed lines are negative. 







Fig. 11. — Slices through {t^^ (t^) and (tj'^ at a pivot scale of fc* = 4.4 x 10 ^kdrf) and with «b - 0. The scale on 
{t^"^ has been truncated to show structure and dashed areas are negative. 
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Fig. 12. — Slices through {r^J, (^^y and yr^j at a pivot scale of A;* = 4.4 x 10^^ kdrj) and with hb = -5/2. The scales 
have been truncated to highlight the structure, dashed areas are negative and the inset focuses on the vicinity of the 
degenerate line for {t^Y 



